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ABSTRACT 

We construct a simple model for radioisotopic enrichment of the protosolar 
nebula by injection from a nearby supernova, based on the inverse square law 
for ejecta dispersion. This idealized model concisely and explicitly connects the 
observational data to the key astrophysical parameters; as such, it provides a 
useful tool for parameter studies that can inform the computationally expen- 
sive numerical (magneto)hydrodynamic studies essential for a full solution to 
this problem. We find that the presolar radioisotopes abundances (i.e., in so- 
lar masses) demand a nearby supernova: its distance D can be no larger than 
D/Rss < 66 times the size i?ss of the protosolar nebula, at a 90% confidence 
level, assuming 1 M Q of protosolar material. The relevant size of the nebula 
depends on its state of evolution at the time of radioactivity injection. In one 
scenario, a collection of low-mass stars, including our sun, formed in a group or 
cluster with an intermediate- to high-mass star that ended its life as a supernova 
while our sun was still a protostar, a starless core, or perhaps a diffuse cloud. 
Using recent observations of protostars to estimate the size of the protosolar 
nebula constrains the distance of the supernova to D ~ 0.02 pc to 1.6 pc. The 
supernova distance limit is consistent with the scales of low-mass stars formation 
around one or more massive stars, but it is closer than expected were the sun 
formed in an isolated, solitary state. Consequently, if any presolar radioactivities 
originated via supernova injection, we must conclude that our sun was a member 
of such a group or cluster that has since dispersed, and thus that solar system 
formation should be understood in this context. In addition, we show that the 
timescale from explosion to the creation of small bodies was on the order of 1.8 
Myr (formal 90% confidence range of to 2.2 Myr), and thus the temporal chore- 
ography from supernova ejecta to meteorites is important. Finally, we can not 
distinguish between progenitor masses from 15 to 25 M in the nucleosynthesis 
models; however, the 20 M model is somewhat preferred. 

Subject headings: Sun: abundances; Sun: general; stars: formation; stars: cir- 
cumstellar matter; stars: pre-main sequence; nuclear reactions, nucleosynthesis, 
abundances 
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1. Introduction 

The picture of star formation that has been developed over the last few decades provides 
a good broad-brush explanation of isolated, low-mass, Sun-like star formation (e.g. Shu et al. 
1993). However, is that picture valid for the formation of the Sun? The scenario is probably 
unlikely as there is good evidence that the majority of Sun-like stars form in clusters (e.g., 
Carpenter 2000; Lada & Lada 2003). In addition, the theoretical view that protostars form 
in relative isolation from their molecular cloud neighbors is also weakening (e.g., Clark et al. 
2005). 

To better explain the origin of the Sun and Planets, it is important to first understand 
the initial conditions in which the Sun formed. Did the Sun form in isolation or in a small 
group or cluster of stars that has since dispersed? The concept of a lonely Solar System 
birth has been challenged by the idea that the Sun was born near at least one massive star, 
in a group or cluster of stars, based on the study of short-lived radioisotopes in meteorites 
(e.g. Hester et al. 2004). 

The main evidence for short-lived radioactive species is from isotopic anomalies of daugh- 
ter species in primitive meteorites (for recent reviews, see Goswami & Vanhala 2000; Meyer 
& Clayton 2000; Wadhwa et al. 2006). The isotopes are now extinct, but they must have 
been live, and in relatively high abundances, during the early stages of planet formation in 
our Solar System. Possible origins for the measured isotopic abundances are that the Sun 
formed in a location that was enriched either by a stellar wind or a recent core-collapse su- 
pernova (e.g. Lee et al. 1977), hence in or near a cluster or stellar group environment. Other 
possibilities for the short-lived radioactive isotopes include irradiation in the circumstellar 
disk (e.g., Lee et al. 1998) or mixing of ejecta from core-collapse supernova in the interstellar 
medium (ISM) (e.g., Meyer & Clayton 2000; Hester & Desch 2005; Ouellette et al. 2005). 
However, new results reveal that the early Solar System contained significant amounts of the 
short-lived 60 Fe radionuclide (e.g. Tachibana & Huss 2003; Mostefaoui et al. 2005) that can 
not be explained with the irradiation models or ISM mixing models (e.g. Tachibana & Huss 
2003; Gounelle et al. 2006). 

In this paper, we make the assumption that a supernova event occurred near the early 
solar nebula, such that short-lived radioactive species were injected into the solar nebula 
(e.g. Cameron & Truran 1977; Goswami & Vanhala 2000). Although, the 60 Fe abundance 
could also be explained by ejecta from Asymptotic Giant Branch (AGB) stars, novae, or 
Type I supernovae, the likelihood of such late stellar stage objects overlapping with a star 
formation region is much smaller than a supernova that was born in the same cluster (e.g. 
Kastner & Myers 1994; Hester & Desch 2005). 
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With the assumption of a nearby supernova, we construct a simple model that directly 
and explicitly relates the observed pre-solar radioisotope abundances to the supernova dis- 
tance and the size of the protosolar nebula into which they were injected. Our model, while 
idealized, clearly illustrates that the radioisotope abundances directly constrain (and in par- 
ticular, place an upper limit on) the ratio of the supernova distance to the size of the nebula. 
This type of model has been discussed before (e.g. Vanhala & Boss 2000, 2002; Ouellette 
et al. 2005) with more emphasis on the detailed hydrodynamics of the injection of supernova 
material into the nebula. In this paper, we focus on a statistical parameter study, with 
emphasis on the supernova distance D. Our simple model allows for a rapid search of the 
parameter space and shows that in most scenarios, a nearby supernova (i.e., D < a few pc, 
and <C 1 kpc) is required if any radioisotopes came from a supernova. This has important 
qualitative implications for the nature of the Sun's formation, and quantitatively our results 
can inform and calibrate the more computationally expensive hydrodynamics calculations 
needed to solve the problem in all of its detail. 



2. Supernova Radioisotope Injection into the Protosolar Nebula 

The abundance of radioisotopes injected by a supernova explosion depends on-and thus 
probes-the properties of the Type II core-collapse supernova and its explosion, the propaga- 
tion of the supernova remnant and entrained ejecta, and the ability of the protosolar nebula 
to capture the ejecta that impinge on it. Each of these aspects of the problem are known to 
be quite complex in its details (e.g. Chevalier 2000; Meyer & Clayton 2000; Vanhala & Boss 
2000, 2002), and a full understanding of the problem must embrace this richness. Nonethe- 
less, a simplified analytic treatment of the parameter space of the problem not only offers 
physical insight, but also provides a benchmark for comparison with the more elaborate 
descriptions. 



2.1. Radioactivity Distance 

In this spirit, we adopt the following idealized model (also see, Vanhala & Boss 2000; 
Ouellette et al. 2005). We assume a supernova event occurred at a distance D from the 
early solar nebula. The ejecta from the supernova blast is then spread with a "flux" (in fact, 
surface density) that follows the inverse square law. It then follows that the mass fraction 
Xi = Mgs,i/Mss, of a radionuclide i, observed in the Solar System is 

i_/mj M ss 47rL> 2 1 j 
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where M S s,i is the mass of radionuclide i in the early solar system, M ss is the effective "target" 
or "absorbing" mass of the solar nebula, M SN>i is the amount of radionuclide i produced in 
the supernova, and i?ss is the solar system radius. The exponential term accounts for the 
decay of the radionuclide, with decay constant Tj, in the time interval t between the supernova 
explosion and solar system formation (more precisely, until meteorite formation). Finally, the 
/inj < 1 measures the "injection efficiency," i.e., the ratio of incident supernova radioisotope 
debris that is actually incorporated into the nebula. For example, Vanhala & Boss (2000) 
discuss the injection process in detail suggesting that in some cases the injection occurs via 
Rayleigh- Taylor instabilities with an efficiency of approximately / ~ 0.1. 

This model is simplified, and thus idealized, by assuming not only spherical, homogenous 
ejecta from isotropical supernova ejecta, but also that the supernova ejecta will be completely 
integrated and well mixed in the solar system material. On the other hand, the model readily 
provides insight into the proximity of the supernova that ejected the short-lived isotopes into 
the early solar nebula. 

From Equation 1, we can derive a measure of "radioactivity distance", which is the ratio 
of the supernova distance to the solar nebula radius. 1 



where the fiducial values are appropriate for 60 Fe and 1 M of solar nebula material. In 
other words, the radioactivity distance is only dependent on four variables. Arguably, the 
two with the largest uncertainties are the supernova yields MsN,i and presolar radioisotope 
abundance Xj. The supernova yields Msn,« are gathered from models of nucleosynthesis in 
massive stars (e.g., Woosley & Weaver 1995; Rauscher et al. 2002), which depend chiefly on 
the mass of the progenitor star. It is interesting to note that the input parameters appear 
under a square root in the expression for radioactivity distance, which is correspondingly 
less sensitive to errors in the inputs. 

The abundance at isotopic closure of a once-live radionuclide parent % V is obtained 
from meteoritic measurement of anomalous abundances in its daughter product l V, where 
the decay scheme is % V — > l T> (see recent review, Wadhwa et al. 2006). This procedure yields 
the initial ratio l V/V ie f of the radioactive parent nuclide to a stable (and abundant) reference 
isotope V re { of the same atomic species. For example, studies of 60 Ni excesses in meteorites 



lr The radioactivity distance is a close analog of a luminosity distance: radioisotope nucleosynthetic yield 
MsN,i takes the role of luminosity and solar abundance (in fact, surface density Mss.i/tt-Rss) takes the role 
of observed flux. 




1/2 
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give a measure of the pre-solar 60 Fe/ 56 Fe ratio (e.g. Tachibana & Huss 2003; Mostefaoui et al. 
2005; Tachibana et al. 2006). 

Table 1 lists the short-lived radionuclides in the early solar nebula that are most likely to 
be contaminated by a Type II supernova ejecta (e.g. Wadhwa et al. 2006). We have excluded 
92 Nb and 244 Pu since the former has very large uncertainties and the latter is not predicted 
in the two models used in our analysis. The ratios used are for the initial solar system 
abundance, which is defined at the creation of the calcium-aluminum-rich inclusions (CAIs), 
not necessarily at isotopic closure. This, however, creates difficulties when referencing the 
age at isotopic closure to the initial solar system abundance for some radionuclides. The 
36 C1/ 35 C1 and 41 Ca/ 40 Ca ratios are impacted the most from this uncertainty and may be 
considered lower limits (e.g., Wadhwa et al. 2006). 

To recover the observed mass fraction for the presolar nebula (Table 1), we use the 
expression 

Xi = A * X h (3) 

where Ai is the atomic weight of the isotope. The early solar V TC {/B. abundance and the 
hydrogen mass fraction X H = 0.711 were taken from Lodders (2003), Table 2, where the 
associated uncertainties from the meteoritic ratios and abundances have been propagated. 
In most cases, the largest uncertainty is from the meteoritic ratios; the exception of which 
is the radionuclide ratio 129 I/ 127 I, which has a large abundance uncertainty. 

For the supernova radioisotope yields Msn,i, we use the results of theoretical models. 
To illustrate the (significant) uncertainties in these difficult calculations, we will compute 
results using both the model of Woosley & Weaver (1995) (at solar metallicity), which spans 
the mass progenitor range 11 — 40 M , and that of Rauscher et al. (2002), which updates 
much of the physics of the earlier model, but spans a narrower mass range. The model 
uncertainties are clearly large: for a given progenitor mass, the yields for the same isotope 
in the two models can diverge in some cases by as much as a factor of ten. We thus believe 
it is not overly conservative to adopt an uncertainty on the model yields as o"M SN i = -^SN,b 
i.e., a factor of two. 

In addition, although a supernova can probably occur with stellar masses of 8 to 10 
Mq, and indeed the lower mass supernova may be the most likely for the case of the Sun, 
the fate of these stars, and their nucleosynthesis, remains unclear. The main reason is 
that the lower mass objects are likely to behave differently than the standard core collapse 
SNe with more associated r-processed material (Wheeler et al. 1998). Unfortunately, these 
systems are not well studied, so not further discussed in this paper. Only the mass range 
of 11 — 25M is considered, which provides results for low- mass supernovae and an overlap 
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of the two nucleosynthesis models. The 30, 35, and 40 M models in Woosley & Weaver 
(1995) provide flat yields in comparison to the 25 M & model, with the exception of 53 Mn 
that varies significantly depending on the assumed kinetic energy. 

2.2. Results 

In this section, we present numerical results for the radioactivity distance ratio D/R ss . 
We will assume a perfect injection efficiency /i n j = 1 and a solar system nebula mass Mss = 
1M Q , in both the text and figures, and we will revisit these assumptions in our discussion. 

1 /2 

Note that the results can be scaled for other efficiencies or masses via D/Rss oc / in j or 
D/Rss oc Mss' 1 ^ 2 - In §3, these results will be applied to two basic cases: (1) injection into 
a core or circumstellar envelope, where f\ n j ~ 0.1, Mss — 1-W©, and -Rss — 0.05 pc; and (2) 
injection into a circumstellar disk, where /i n j ~ 1, Mss = Mdisk — O.O1M , and Rss — 100 
AU. 

Figure 1 shows the radioactivity distance for the nine short-lived radioactive isotopes 
tabulated in Table 1, for the case t/Vj <C 1 of negligible delay between radioisotope nucle- 
osynthesis and incorporation into meteorites. Both nucleosynthesis models are shown: the 
Rauscher et al. (2002) model with solid lines and the Woosley & Weaver (1995) model with 
dashed lines. As the figure is already complicated, the error bars on the Woosley & Weaver 
(1995) model are not drawn. Broadly speaking, the order of magnitude of the results are 
in line with the fiducial estimate in Equation (2), spanning D/Rss ~ 10 1 — 10 3 . We im- 
mediately see that if any pre-solar radioisotopes originated in a supernova, then (1) the 
explosion had to be close, within 10-10000 solar nebula radii; and (2) the supernova had to 
be recent, t/Vj < 1, else many of the shortest lived radioisotopes would have decayed and/or 
the required distance would become unfeasibly close. 

In more detail, the large range of radioactivity distance in Figure 1 illustrates the com- 
plexity of the problem and the uncertainties in such a model. The nucleosynthesis uncer- 
tainties alone are clearly very large, as the two supernova models lead to distance estimates 
that can vary by factors as high as ~ 3. In the ideal case, in which all radioisotopes did 
originate in a supernova as described in the model, one would expect the curves for all iso- 
topes to converge to a single unique value at some progenitor mass; this is not the case. In 
fact, it is clear from Figure 1 that the mass of the progenitor is unconstrained. Nonetheless, 
the majority of the isotopes are clustered near the middle of the graph, with some isotopes 
preferring larger radioactive distances ( 41 Ca and 53 Mn) and smaller radioactive distances 
( 146 Sm). 
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Two possible explanations for the stratification in Figure 1 suggest themselves: mul- 
tiple sources of short-lived radioisotopes and decay of the radioisotopes. We have thus far 
considered the case of a supernova as the unique source of presolar radioactivities, but other 
sources certainly are possible and indeed must be present at some level. One source is the 
abundance of radioisotopes in the general interstellar medium, which achieves a steady state 
between production and decay. The steady-state abundance is essentially the production 
abundance, reduced by the ratio Tj/r* of decay time to a star-formation time r* appropri- 
ate for the global ISM. Consequently, the highest steady-state abundances will be those of 
long-lived species, such as 53 Mn, 146 Sm, 107 Pd, and 129 I and others (Meyer & Clayton 2000). 

Another more promising source of presolar radioisotopes is in situ production through 
protostellar cosmic rays in the circumstellar disk ( 41 Ca and 53 Mn; Lee et al. 1998). Indeed, the 
presence of presolar 10 Be (McKeegan et al. 2000) suggests that this may well be the case, since 
this isotope is in general only produced by cosmic rays, and not at all by supernovae (e.g., 
Vangioni-Flam et al. 1996). We note that there should also be a steady-state 10 Be abundance 
in the general ISM. If cosmic-ray bombardment occurs over a timescale 3> t\q = 2.2 Myr, 
then immediately afterwards, the 10 Be abundance relative to the stable isotope is 10 Be/ 9 Be ~ 
r io/ r * ~ 10~ 3 (1 Gyr/r*). To explain the full pre-solar value 10 Be/ 9 Be ~ 10~ 3 (McKeegan 
et al. 2000) would seem to require an effective Galactic cosmic-ray nucleosynthesis timescale 
of ~ 1 Gyr, which is short compared to the ~ 9 Gyr cosmic age at the solar birth. 
In addition to this Galaxy-wide beryllium production, the nascent Solar System was likely 
subject to transient and localized enhancements in Galactic (i.e., high-energy) cosmic rays 
that would lead to some level of enhancement in spallation radioisotope production (Fields 
et al. 1996). We note that a nearby supernova represents a local cosmic-ray accelerator. In 
addition, as the protosolar core collapses, it becomes opaque to the already high flux of 10 Be 
Galactic cosmic-rays and traps them (Desch et al. 2004). Nonetheless, 60 Fe is still the one 
isotope not produced at the needed levels in the cosmic-ray scenario, and thus points to the 
need for at least some supernova material (e.g., Hester & Desch 2005). 

In any case, local production (and/or other extra sources, including the possibility of 
multiple supernovae) would mean that the presolar radioisotope abundance would have 
two, or more, components. In those conditions, the radioactivity distances estimated in this 
paper will, by definition, be lower limits to the supernova distance. To disentangle these 
sources is not a trivial task, but not hopeless either, as we are aided by (1) the wide range 
of radioisotope lifetimes, and (2) the rather different abundance patterns arising from the 
various nucleosynthesis mechanisms. For example, if there is a significant local cosmic-ray 
production of, say, 26 Al, this would reveal itself as an anomalously low radioactivity distance 
for this isotope relative to a supernova-only species such as 60 Fe. Such an analysis is beyond 
the scope of this paper. 
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Based on probable secondary sources, we exclude 146 Sm from the parameter study in 
§2.3. The isotope has the longest half-life in Table 1 by almost a factor of three. As 
suggested by Meyer & Clayton (2000), meteoritic measurements of 146 Sm are most likely 
contaminated by a steady-state component in the ISM maintained by Type la supernovae. 
If the abundance is increased, the derived radioactive distance decreases in Equation 2, which 
explains its location at the bottom of Figure 1. On the other hand, 53 Mn, although within 
the error bars of 60 Fe, demands a high radioactive distance. Yet, like 146 Sm, the measured 
53 Mn abundance is expected to have been contaminated by the steady-state ISM abundance 
(Meyer & Clayton 2000), implying an even larger distance for the supernova. However, 53 Mn 
is produced during explosive silicon burning, which occurs at the very end of the supernova 
process. Because the Type-II supernova explosion mechanism is still poorly understood, it 
is difficult to accurately predict the products of explosive silicon burning (Woosley et al. 
2002). In addition, the overabundance of 53 Mn in the supernova ejecta has been explained 
by the probable fallback of the innermost layers of the star (Meyer 2005). For these reasons, 
we also exclude 53 Mn from the parameter study in §2.3. 

Finally, within the context of the model in this paper, the radioactivity distance es- 
timates for each species will vary due to the effect of radioactive decay, if there is any 
significant time delay between nucleosynthesis and incorporation into meteorites. Of course, 
for any material injected by a nearby supernova, some time must inevitably elapse between 
radioisotope production and incorporation into meteorites, in particular the CAIs. At mini- 
mum, there must be a nonzero time-of-flight between the supernova and the presolar nebula; 
also, the time for creation of small bodies should be ~ 1 to 10 Myrs based on circumstellar 
disk lifetimes (e.g., Hartmann 2005). If the supernova impact itself initiated the collapse or 
merely seeded a nearby cloud of the presolar nebula (Cameron & Truran 1977), then there 
is an additional delay due to the collapse, which is estimated to take ~ 0.1 to 10 Myrs (e.g., 
Shu et al. 1993; Vanhala & Cameron 1998; Vanhala & Boss 2000; Tassis & Mouschovias 
2004). 

Any delay would imply that the supernova radioisotope signal is underestimated, and 
correcting for this effect decreases the radioactivity distance in Equation (2). Thus, the 
un-delayed results of Figure 1 represent upper limits to the distance. However, adjustment 
to the radioactive distance is exponentially sensitive to the isotope lifetime; the shorter-lived 
species such as 41 Ca and 36 CI will vary strongly for modest delays of ~ 100 kyr, while the 
longer-lived species (r > 5 Myr) will be essentially unaffected over timescales appropriate 
for the full ensemble of isotopes. For example, this explains why 41 Ca demands the largest 
radioactivity distance in Figure 1; 41 Ca has decayed the most. On the other hand, because of 
their short half-lives the initial 41 Ca and 36 CI abundances are also among the most difficult 
to determine (e.g., Wadhwa et al. 2006). Indeed, their estimated abundances in Table 1 are 
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really lower limits, which implies upper limits in Figure 1. 



2.3. Parameter Modeling 



In order to comment on the degree of concordance in Figure 1, we define the x 2 01 the 
data at each progenitor mass with assumed radioactivity distance (D/R ss ) gucss to be 



We evaluated a grid of x 2 values for (D/R ss ) guC ss ranging from 1 to 200, in steps of I. For 
each nucleosynthesis model and progenitor mass, a separate x 2 value was calculated. In 
addition, a decay time ranging from to 2.5 Myrs, in steps of 0.1 Myrs, was included. This 
gives a total of 5200 possible combinations for each progenitor mass. A parameter set was 
considered an acceptable fit at a 90% confidence level, i.e. likelihoods > 10%. 

Due to the large uncertainty in the data and nucleosynthesis models, there is a large 
range of acceptable fits. The Rauscher et al. (2002) supernova nucleosynthesis model, with 
the larger number of predicted isotopes but smaller mass range, has the largest constraints on 
the parameters. Therefore, we will discuss those fits in detail. Figure 2 shows the calculated 
likelihoods of the parameter study. The lowest contours in Figure 2 are 10%. 

With accepted likelihoods greater than 10%, the mass is not well constrained (i.e. all 
mass models fit); however, it is clear that the mass model of 20 M is the highly preferred 
model. In fact, the best fit (likelihood of 93%) is with a progenitor mass of 20 M , decay time 
of 1.8 Myrs, and a D/R55 ratio of 22. However, it is important to keep in mind that all fits 
in Figure 2 are statistically equivalent in our 90% confidence level. Although the preferred 
X 2 of the 20 M model is perhaps intriguing and should be explored more in detail, it is 
probably a conspiracy of terms. The tangible reason that the 20 M Q model fits the data 
so well is due to the increased predicted abundance of 41 Ca and 36 CI in that model. It is 
crucial to remember that the strict meteoritic data for 41 Ca and 36 C1 are upper limits, which 
can heavily impact the result by decreasing the decay time necessary to achieve an overlap. 
Figure 3 shows the radioactive distance at four different decay times, t = 0, 0.1, 1.8, and 2 
Myrs, t = 1.8 Myrs was included as it is the best fit time. The figure illustrates the evolution 
with different decay times. 

It is clear from the above discussion that the most constrained aspect of the simplified 
problem is the radioactivity distance with some constraint on the decay time. Using a 90% 




(4) 
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confidence level, the distance ratio and decay time is constrained to be 



5 < D/Rss < 66 



(5) 



and 



< t < 2.2 Myrs, 



(6) 



respectively, using the Rauscher et al. (2002) nucleosynthesis model. In order to compare 
the distance ratio to other distance estimates, we have to use assumptions for the early solar 
system mass and radius. The large time range is unsurprisingly consistent with many other 
estimated time scales (e.g., Meyer et al. 2003; Meyer 2005; Ouellette et al. 2005). 



The largest assumption in Figure 3 is that a supernova exploded nearby the protosolar 
nebula such that up to 2.2 Myrs later those isotopes were incorporated into meteorites. The 
simplest and most probable explanation is that our Sun formed in a loose cluster containing 
a massive star that later dispersed. In that case, was the Sun triggered by the supernova 
(e.g., Cameron & Truran 1977), or did it form at a similar time as the progenitor star (e.g., 
Chevalier 2000)? 

One must keep in mind that the temporal choreography is crucial to this discussion. 
Our results show that the time scale for the isotopical inclusion into solid objects ranges 
from to 2.2 Myrs, including travel, grain growth, and planetesimal formation. That is a 
particularly quick timescale when compared to the lifetime of a protostar. In addition, if 
we consider that the expected timescale for large body formation occurs at the end of the 
protostellar lifetime, during the Class II/III stage or later (e.g., Hartmann 2005), then the 
likelihood of this particular supernova triggering our solar system does exist, but it must 
evolve quickly. Nonetheless, we discuss these scenarios using the current best estimates of 
early nebula radius and mass. 



In the case of siblings, the protosun and the progenitor star are formed on similar time 
scales (within ~ 1 Myrs, Lada & Lada 2003): possibly a birth case similar to the cluster of 
about 1600 low-mass young stars in the Orion Nebula within ~2 pc of the the Trapezium 
group (core radius of 0.2 pc, Hillenbrand & Hartmann 1998). Within the inner 0.4 pc, 85% 
of the stars have L-band excesses that suggest circumstellar disks (Lada et al. 2000). In 



3. Implications for the Presolar Environment 



3.1. Siblings at Birth 
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spite of the harsh environment near the massive cluster, the low-mass young stars near the 
central star 8 1 Ori C have observable disks due to the photoionization effects of the massive 
star (e.g., McCaughrean & O'dell 1996). 

Clearly from Figures 2 and 3, the isotopic evidence does not well constrain the mass of 
the progenitor star. Nonetheless, the main sequence lifetime of the progenitor star can be 
estimated from recent models (e.g., Romano et al. 2005) as 7 to 10 Myrs for the mass range 
in Figure 3. In addition to the main sequence lifetime, one must include the protostellar 
lifetime of the massive star. Although this value is poorly known, the best estimates range 
from 10 4 to 10 5 years (e.g., Yorke 1986; McKee & Tan 2002), a fraction of the main sequence 
lifetime. 

Although the time frame of ~10 Myrs is significant, low-mass stars evolve at a much 
slower pace than higher mass stars. In fact, the longest stage may be before collapse even 
begins on the order of the ambipolar diffusion timescale, which ranges from 0.1 to 10 Myrs, 
depending on the initial mass-to- flux ratio and degree of ionization in the cloud (e.g., Tassis 
& Mouschovias 2004). That in combination with the standard spread of ages in clusters of 
~1 Myrs (Lada & Lada 2003) implies that the protostellar nebula could have been in many 
possible evolutionary stages (also see Hester & Desch 2005): (1) the youngest protostellar 
evolution stage, the so-called Class stage, which is dominated by an circumstellar envelope 
that implies .Rss ^ 5000 AU and M ss ~ 1 M (e.g., Looney et al. 2000, 2003). Although 
the injection efficiency is probably closer to 10% (Vanhala & Boss 2002), we conservatively 
use /i n j ~ 1 (i.e., larger distance); or (2) an older protostar, Class I/II, dominated by the 
circumstellar disk that implies Rss — 100 AU and Mss — 0.02 M (e.g., Looney et al. 2000), 
which is similar to the minimum mass solar circumstellar disk (e.g., Weidenschilling 1977). 
We conservatively assume a large injection efficency of / in j ~ 1. 

From our analysis in §2.3, the radioactivity distance ranges from 5 to 66. Using that 
range of radioactivity distance ratio fits with the above size and mass estimates, yields 
distances of 0.12 to 1.6 pc and 0.02 to 0.22 pc, respectively (cf. Chevalier 2000; Ouellette 
et al. 2005). Note that a very low disk mass in older Class III sources or even debris disks 
would result in a supernova distance that is many pc away in this formalism (i.e. M ss is 
very small). However, it is unclear how those systems could manufacture the meteorites 
necessary for the short lived isotope detection. The ejecta must be incorporated into the 
objects at an early enough time to be detected in meteorites. It is also important to note 
that these derived distances are very consistent with low-mass systems seen in clusters with 
large mass members (e.g., Orion or Eagle Nebula). 
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3.2. Triggered Birth 

If we assume that a supernova event injected short lived isotopes into the early solar 
nebula, one must also consider that this event could have triggered the collapse of the nebula 
(e.g., Cameron & Truran 1977). Many recent models have shown that moderately slow shock 
waves (20-45 km/s) can trigger the collapse of cores (e.g., Boss 1995; Vanhala & Cameron 
1998): too little and no collapse, too much and the core is torn apart. A quick collapse process 
(~ 10 4 — 10 5 yrs, Vanhala & Cameron 1998) may also account for the low radioactive distance 
(i.e. large Xj) of the majority of the longer half-life isotopes ( 146 Sm, 129 I, 182 Hf, and 107 Pd) 
seen in Figure 1; the molecular cloud core in which the Sun formed would not have time to 
decay the initial ISM component. On the other hand, Hester & Desch (2005) argue that a 
supernova triggering a collapse and injecting that collapse with its radioisotopes, is probably 
less common than the scenario where supernova ejecta intercept a protostar. 

In the triggering scenario, the solar nebula would not yet have evolved to a Class 
protostar, implying two possible states: (1) a starless core, which has Rss — 7000 AU and 
typical mass of Mss — 2 M (e.g., Kirk et al. 2005). Again, although the injection efficency 
is probably closer to 10% (Vanhala & Boss 2002), we conservatively use /i n j ~ 1 (i.e., larger 
distance); or (2) a quickly triggered (e.g. < 1 Myrs), diffuse cloud as modeled by a Bonner- 
Ebert sphere, considered by Vanhala & Boss (2000, 2002), with R ss ~ 0.06 pc, mass of 
M ss ~ 1M Q , and / in j ~ 0.1. With these two possibilities, the estimated ranges of distances 
are 0.12 to 1.6 pc and 0.06 to 1.2 pc (or 0.3 to 4 pc for the overly conservative case of / in j ~ 1), 
respectively. This can be compared to the values of Vanhala & Boss (2000, 2002), 5.2 to 26 
pc. Vanhala & Boss (2000, 2002) comment on the mismatch of distance and offer possible 
resolutions. On the other hand, the possibility of a supernova inducing a quick collapse of a 
very diffuse cloud, i.e., Rss > 0.1 pc, will significantly increase the derived distance to the 
supernova. 



3.3. Travel Time 

With an assumed ejecta speed of 30 km/s and the maximum estimated distances from 
above (e.g., 1.6 pc or even most conservatively 4 pc), the travel time can be estimated at 
~ 5 x 10 4 and 1 x 10 5 yrs, respectively. Using our best fit decay time of 1.8 Myrs, this 
implies that the ejecta material mixed with the solar nebula material until meteorites were 
produced. Given nearly 1.8 Myrs of residency within the solar nebula, it is possible that 
the SN ejecta would become homogenized in the cloud. This would be consistent with the 
formation timescale of calcium-aluminum-rich inclusions (CAIs), on the order of a million 
years (e.g., Wadhwa & Russell 2000). 
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4. Conclusions 

We present a simple model relating the pre-solar abundances of short-lived radioiso- 
topes to the properties of the injecting supernova event. While our model is idealized, it is 
also explicit and computationally simple, yielding a parameter study whose general aspects 
are intriguing. In general, we are driven to conclude that a very nearby supernova event 
must be invoked if any presolar radioactivities are due to injection from a single explosion. 
Specifically, to a formal confidence level of 90%, the distance from the Solar Nebula and the 
supernova was in the range of 0.1 to 1.6 pc, for siblings and up to approximately 4 pc (in 
the most conservative, and unlikely, case of perfect injection) for triggered birth, using the 
Rauscher et al. (2002) nucleosynthesis models for 15 to 25 M progenitors. This sounds sur- 
prisingly close, but it is consistent with typical distances found for low-mass stars clustering 
around one or more massive stars. We posit that our Sun was a member of such a cluster 
that has since dispersed. 

It is important to note that our distance estimates imply that the Sun was directly 
influenced by the massive star's photoionization effects, much like proplyds in Orion. To 
understand the Sun's formation, one must include the large role played by such a complicated 
environment. The minimum distance is also interesting, as it is similar to the distance at 
which the disk would probably be destroyed < 0.25 pc (Chevalier 2000). Our estimated 
distance range includes the expected uncertainties of the nucleosynthesis models and the 
radioisotopic measurements. Although the physical parameters of the solar nebula also play 
a large role in these estimates, we use nearly all protostellar evolutionary states to compile 
our range. 

In addition to the distance of the supernova, we constrain the time scale of the explosion 
to the creation of small bodies to a 90% confidence range of to 2.2 Myrs. The majority of 
this time was probably spent in our solar system during the building of small objects from 
the circumstellar dust; in fact, it is consistent with the estimates for chondrule production. 
It is important to note that the temporal choreography from supernova ejecta to meteorites 
is crucial. The upper limits of our fit decay times make it possible for the supernova to have 
triggered the fast collapse of our Sun, seeded it with short-lived isotopes, and created CAIs, 
but it seems more likely that the Sun was already a protostar during the explosion (e.g., 
Hester & Desch 2005). 

In addition to the inherent uncertainties in our model, there are additional assumptions 
that would imply a smaller or larger supernova distance. There are two clear cases that 
would cause our deduced distances to be an upper limit: (1) an inclined circumstellar disk 
(Class I/II source) would have a lower cross-sectional area and not collect the full supernova 
contaminates as a face-on disk, and (2) the ejecta may not be effectively stopped by the 
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circumstellar disk as assumed (/ in j <C 1) (i.e., effectiveness of the injection process; see 
Vanhala & Boss 2002). Similarly, there are two opposing cases that would make our deduced 
distances lower limits: (1) highly inhomogeneous ejecta (i.e. bullets as seen in Cas A, 
Willingale et al. 2003) would not allow a predictable distance, and (2) although the material 
is injected uniformly onto the circumstellar disk, the inner regions of the disk has a higher 
density than the outer regions (e.g., Mundy et al. 1996), so if there is no mixing, the measured 
meteorites are produced in situ with abundances lower than expected. These crucial issues 
lie within the domain of sophisticated (magneto) hydrodynamical simulations that will yield 
a fuller understanding of this problem; our simple model serves as a tool that highlights and 
quantifies the issues that the full numerical models address. 

In closing, we note that the analysis presented here for extinct presolar meteoritic ra- 
dioactivities closely parallels considerations of live geological radioactivities as signatures of 
relatively recent nearby supernovae (e.g., Ellis et al. 1996; Knie et al. 2004). Both analyses 
leverage the uniqueness of radioisotopes as a signature of supernova activity, and both use 
the inverse square law in the same way to gauge distance. Moreover, Fields et al. (2005) 
emphasize that given a confirmed nearby supernova event, one can then turn the problem 
around, reading the presence or absence of, and detailed abundances of, all radioisotopes 
as direct probes of supernova ejecta and thus nucleosynthesis. This adds another way that 
meteoritic anomalies serve as fossil testimony to our massive sibling. 
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NASA under Origins grant No. NN06GE41G. B.D.F. acknowledges support by the National 
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Table 1. Adopted Short-lived Radioactive Isotopes 



Radio- a 


Ref. 


*l/2 


Meteoritic b 


Meteoritic b 


Reference 


Reference 


Mass 




Mass 




isotope 




Ratio 


Uncert. 


Abundance 


Uncert. 


Fraction 




Uncert. 




i-p 


V re{ 


(Myr) 


l v/v iet 


ff.p 


Pref/H 




Xi 








26 Al 


27 A1 


0.72 


5.9 x 10" 5 


0.6 x 10" 5 


3.46 x 10" 6 


0.16 x 10" 6 


3.77 x 10- 


9 


0.42 x 10" 


-9 


36 CI 


35 C1 


0.30 


1.6 x 10" 4 


0.7 x 10" 4 


2.15 x 10" 7 


0.29 x 10" 7 


8.82 x 10" 


10 


4.04 x 10" 


-10 


41 Ca 


40 Ca 


0.10 


1.4 x 10" 8 


0.1 x 10" 8 


2.59 x 10" 6 


0.18 x 10" 6 


1.06 x 10" 


12 


0.10 x 10" 


12 


53 Mn 


55 Mn 


3.7 


2.8 x 10" 5 


0.3 x 10" 5 


3.77 x 10" 7 


0.26 x 10" 7 


3.98 x 10" 


10 


0.51 x 10" 


-10 


60 Fe 


56 Fe 


1.5 


7.5 x 10" 7 


2.5 x 10" 7 


3.45 x 10" 5 


0.24 x 10" 5 


1.10 x 10" 


9 


0.38 x 10" 


-9 


107p d 


108p d 


6.5 


2.0 x 10" 5 


0.2 x 10" 5 


5.90 x 10" 11 


0.41 x 10" 11 


8.98 x 10- 


14 


1.09 x 10" 


-14 


129j 


127j 


15.7 


1.1 x 10" 4 


0.1 x 10" 4 


4.10 x 10" 11 


1.13 x 10" 11 


4.14 x 10" 


13 


1.20 x 10" 


-13 


146 Sm 


144 Sm 


103 


7 x 10" 3 


2 x 10" 3 


1.05 x 10" 11 


0.10 x 10" 11 


7.60 x 10" 


12 


2.28 x 10" 


-12 


182 Hf 


180 Hf 


8.9 


1.1 x 10" 4 


0.1 x 10" 4 


6.99 x 10" 12 


0.64 x 10" 12 


9.95 x 10" 


14 


1.28 x 10" 


-14 



a Only those isotopes suspected to come from core-collapse supernovae (e.g., Wadhwa et al. 2006) are included, except 
92 Nb and 244 Pu. 92 Nb and 244 Pu were omitted because of a very high uncertainty and a lack of predictions in the two 
supernovae ejecta models used in this paper, respectively. The ratio given here is the initial solar system abundance, which 
is measured at the creation of the CAIs, not necessarily at isotopic closure. 

b Other references used are 26 Al: Young et al. (2005a,b); Cosarinsky et al. (2006); 36 C1: Lin ct al. (2005) (The ratio used 
is dependent on the 26 A1/ 27 A1 ratio in the measured mantle compared to the CAI. The ratio used was consistent, within 
the error bars, with our assumed value. However, the 26 A1/ 27 A1 ratio in the mantle was not well detected, so this value is 
more correctly taken as a lower limit.); 41 Ca: Srinivasan et al. (1996) (This ratio was measured in a CAI with a measured 
26 A1/ 27 A1 ratio that was consistent, within the error bars, with our assumed value. However, it is important to note that if 
that particular CAI formed late, even by a small amount, the measured 41 Ca abundance ratio would be heavily impacted. 
Again, this ratio is more correctly taken as a lower limit.); 53 Mn: Nyquist et al. (1999); 60 Fc: Tachibana ct al. (2006); 
1C)7 Pd: Hauri et al. (2000); Carlson & Hauri (2001) (There is a suggestion that the initial abundance ratio may have been 
as high as 40 X 10" 5 .); 129 I: Brazzlc (1999); 146 Sm: Lugmair & Galcr (1992); 182 Hf: Klcinc ct al. (2005, 2006). Often the 
published error bars are 2cr, but to be conservative we used the errors as la. If the estimated error was less than 10%, we 
used a conservative 10% error. 

c From Table 2 Lodders (2003), using the initial solar system abundances, i.e. correction for gravitational settling. 
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Progenitor Mass (M Q ) 



Fig. 1. — The radioactive distance ratio without any decay time. The supernova production 
models of Woosley & Weaver (1995), dotted lines, and Rauscher et al. (2002), solid lines, are 
used. The estimated error for each data point is shown only for the Rauscher et al. (2002) 
models in order to minimize confusion in the Figure. 




Fig. 2. — The likelihoods of the x 2 parameter search for 5 of the models in Rauscher et al. 
(2002). The mass of the progenitor in the model is listed in each panel. The x-axis is the 
radioactivity distance ratio: from 1 to 200 in steps of 1. The y-axis is the time delay from 
ejection to incorporation into CAIs: from to 2.5 Myrs in steps of 0.1. The contours of 
liklihood are in 10 to 90% in steps of 10%. 
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Fig. 3. — The Radioactive Distance ratio with a decay time of 0, 0.1, 1.8, and 2 Myrs for 
the Rauscher et al. (2002) radionucleosynthesis model. Note that 146 Sm and 53 Mn have been 
excluded for reasons described in the text. 



